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15 Summary 

16 1. Understanding how landscape features affect functional connectivity among populations 

17 is a cornerstone of landscape genetic analyses. However, parameterization of resistance 

18 surfaces that best describe connectivity is largely a subjective process that explores a 

19 limited parameter space. 

20 2. ResistanceGA is a new R package that utilizes a genetic algorithm to optimize 

21 resistance surfaces based on pairwise genetic distances and either CTRCUITSCAPE 

22 resistance distances or cost distances calculated along least cost paths. Functions in this 

23 package allow for the optimization of both categorical and continuous resistance surfaces, 

24 as well as the simultaneous optimization of multiple resistance surfaces. 

25 3. There is considerable controversy concerning the use of Mantel tests to accurately relate 

26 pairwise genetic distances with resistance distances. Optimization in Resi stanceGA 

27 uses linear mixed effects models with the maximum likelihood population effects 

28 parameterization to determine AICc, which is the fitness function for the genetic 

29 algorithm. 

30 4. Resi stanceGA fills a void in the landscape genetic toolbox, allowing for unbiased 

31 optimization of resistance surfaces and for the simultaneous optimization of multiple 

32 resistance surfaces to create a novel composite resistance surface. 
33 

34 Key-words: CIRCUITSCAPE, circuit theory, cost distance, gene flow, genetic algorithm, 

35 landscape genetics, least cost path, resistance distance, resistance optimization, resistance surface 
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36 Introduction 

37 First coined in 2003, landscape genetics has experienced rapid growth in both the number of 

38 studies and range of analytical methods utilized (Manel et al. 2003; Storfer et at 2010). This 

39 integrative field draws on landscape ecology, spatial statistics, and population genetics to address 

40 a wide range of questions. One of the primary goals of landscape genetic studies has been to 

41 understand how landscape features affect spatial genetic structure (Manel et al. 2003; Storfer et 

42 al. 2007). Rather than simply assessing isolation-by-distance, many landscape genetic studies 

43 quantify the effective distance (i.e. resistance distance) between spatially distinct sample 

44 locations as a function of the landscape matrix (Spear et al. 2005; McRae 2006). In the absence 

45 of direct observation of movement or dispersal across the landscape, resistance distances are 

46 often interpreted as functional connectivity (e.g., Cushman et al. 2006). Critically, to infer 

47 functional connectivity and resistance distance, an appropriate resistance surface must be 

48 parameterized. As defined by Spear et al. (2010), a resistance surface is a spatial layer that 

49 assigns a value to each landscape or environmental feature, with values representing the extent to 

50 which that feature impedes or facilitates connectivity for an organism. 

51 Resistance values of surfaces have been determined using a variety of methods, 

52 including: habitat suitability models (e.g., Wang et al. 2008), telemetry (e.g., Driezen et al. 

53 2007), and most commonly expert opinion (Murray et al. 2009). Less often, parameterization is 

54 informed by empirical movement studies (e.g., Stevens et al. 2006) or spatial prediction of 

55 ecological processes (Peterman et al. 2014). All of these are acceptable approaches, but each 

56 comes with its own caveats. Of particular concern is the fact that expert opinion often fails to 

57 accurately describe the biological or ecological process(es) being modeled (Shirk et al. 2010; 

58 Charney 2012). Even when biological or ecological processes are known and explicitly modeled, 
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59 there is no guarantee that these processes will relate meaningfully to the movement of genes 

60 across the landscape (Peterman et at 2014). 

61 Ultimately, the assignment and evaluation of resistance values generally covers only a 

62 limited parameter space and remains a trial and error process to determine the best resistance 

63 parameters. Graves, Beier and Royle (2013) attempted to alleviate the subjectivity of resistance 

64 surface parameterization by using a search algorithm to maximize Mantel r correlation between 

65 inter-individual genetic distance and least cost path distance. While their optimization procedures 

66 recovered the maximum Mantel r when it existed, they found that resistance estimates were often 

67 imprecise and much smaller than simulated resistance values. They also found that response 

68 surfaces were quite flat, making identification of a global optimum difficult. In contrast, 

69 Peterman et al. (2014) found well-defined global optima when using Ricker and monomolecular 

70 data transformations in combination with optimization algorithms. However, the optimization 

71 procedure utilized by Peterman et al. (2014) is limited to continuous resistance surfaces (e.g., 

72 temperature, moisture, percent canopy cover) and requires an inefficient search of all possible 

73 data transformations. 

74 There are numerous challenges to optimizing resistance surfaces based on pairwise 

75 (genetic) distance data. Among these challenges is the fact that resistance surfaces can have a 

76 high dimensionality, as in land use, land cover surfaces. Another issue is that there currently is 

77 no closed-form expression to determine the landscape resistance values that describe pairwise 

78 distances, potentially making optimization intractable with gradient-based algorithms. Finally, 

79 landscape features and environmental gradients do not exist in isolation. Therefore, an ideal 

80 solution to resistance surface optimization will be to simultaneously optimize multiple resistance 

81 surfaces to create a composite resistance surface. The Resi stanceGA package for the R 
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82 programming environment (R Core Team 2014) has been developed to address these issues, and 

83 to fill a void in the landscape genetic toolbox. While the initial impetus for this package was 

84 landscape genetic analyses, but any pairwise measures across the landscape (e.g., movement 

85 rates) could be utilized to optimize resistance surfaces. 
86 

87 Description 

88 Genetic algorithms 

89 Genetic algorithms (GAs) represent a powerful and flexible stochastic optimization framework 

90 for finding solutions to both discrete and continuous optimization problems (Holland 1975). 

91 Inspired by biological principles, genetic algorithms create a population of individuals 

92 (offspring) with traits (parameters to be optimized) encoded on "chromosomes". The genotypes 

93 (parameter combinations) of each individual solve the fitness function, and the fittest individuals 

94 from each generation survive to reproduce (Sivanandam & Deepa 2007). The GA evolution 

95 process is facilitated by exploration and exploitation (Scrucca 2013). Exploration of parameter 

96 space occurs through random generation of new parameter values resulting from mutation, as 

97 well as exchange of genetic information through crossover. Exploitation ultimately reduces 

98 diversity in the population by selecting the fittest individuals each generation. The population 

99 continues to evolve until a sufficient number of generations have passed without an improvement 

100 in fitness (Scrucca 2013). 

101 Resistance optimization 

102 ResistanceGA utilizes the general-purpose genetic algorithm from the GA R package (g a 

103 function; Scrucca 2013). Briefly, the optimization proceeds as follows: 
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104 1) The original raster surface is imported into R. If the surface is continuous, it is rescaled to 

105 range of 0-10, preserving the relative spacing between all levels. 

106 2) The evolution process starts by generating a random initial population of size n. If a 

107 continuous surface, the selected parameters determine (a) which of eight transformations 

108 will be applied (Fig. 1); (b) the shape of the transformation; (c) the maximum value of the 

109 transformation. If a categorical surface, each level of the resistance surface is reclassified 

110 to the values of the parameters. 

111 3) Using the spatial locations where genetic samples have been collected, either 

112 CIRCUITSCAPE (McRae et al. 2008; McRae & Shah 2009) is called from R to calculate 

113 pairwise resistance distances across the landscape created in step 2 above, or cost 

114 distances are calculated from least cost paths using the R package gdi stance (van 

115 Etten 2014). 

116 4) A linear mixed effects model with a maximum likelihood population effects 

117 parameterization (MLPE) is fit to the data. Pairwise genetic distance is the response and 

118 the scaled and centered pairwise resistance distance is the predictor. The MLPE mixed 

119 effects parameterization is used to account for the non-independence among the pairwise 

120 data (Clarke, Rothery & Raybould 2002). 

121 5) The Akaike information criterion (AIC) is obtained from the fitted MLPE model. AICc is 

122 then calculated by penalizing the model for the number of parameters used during 

123 transformation/reclassification. AICc is the fitness measure that the genetic algorithm 

124 seeks to maximize. Because the ga function works through maximization, the sign of 

125 AICc is reversed. 
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126 6) Steps 2-5 are repeated until the specified number of n individuals have been created. The 

127 genetic algorithm then conducts selection on the population, and the individuals with the 

128 best AICc are carried over to the next generation to form the reproducing population. A 

129 new population is then formed through mutation and crossover. 

130 7) Steps 2-6 are repeated until the specified number of generations have passed without 

131 improvement to the AICc. 
132 

133 Continuous surfaces 

134 There are eight transformations that can be applied to continuous surfaces (Fig. 1). Each 

135 transformation relies on either the Ricker or monomolecular functions (Bolker 2008), as well as 

136 rescaling functions to keep values in positive parameter space . The shape and magnitude of each 

137 transformation are controlled by two parameters, which are adjusted during optimization 

138 (Peterman et al. 2014). During optimization, the genetic algorithm searches different 

139 combinations of transformations, scale parameters, and shape parameters. Linear transformations 

140 are not explicitly included, but all of the monomolecular functions become linear as the shape 

141 parameter increases in value. In this way, linear responses can effectively be modeled without 

142 increasing the number of transformations for the genetic algorithm to test. 

143 Categorical surfaces 

144 Categorical or feature surfaces, such as land cover or roads, can also be optimized using 

145 Resi stanceGA. A surface is considered categorical in Resi stanceGA if it contains 15 or fewer 

146 levels. To make this process tractable, it is necessary to hold the value of one level constant 

147 throughout optimization. Because pairwise resistance values are relative, failing to hold one level 
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148 constant can result in multiple equivalent solutions, and the algorithm may fail to reach an 

149 optimal solution (e.g., relative resistance values of 1, 5, and 10 are equivalent to 2, 10, and 20). 

150 Combining resistance surfaces 

151 Resistance surfaces are simultaneously optimized by sequentially modifying each surface to be 

152 optimized. Continuous and categorical surfaces are each modified as described above, then all 

153 modified surfaces are summed together to create a single, composite resistance surface, which is 

154 then used to calculate pairwise resistance distances or cost distances. 
155 

156 Overview of Resi stanceGA functions 

157 The optimization functions passed to ga rely heavily on the R package raste r (Hijmans 2014) 

158 to import, export and manipulate spatial raster (.asc) files, as well as 1 me4 (Bates et al. 2014) to 

159 fit mixed effects models. To visualize continuous surface transformations, Resi stanceGA uses 

160 the graphing features of ggpl ot2 (Wickham 2009). The functions available in Resi stanceGA 

161 are summarized in Table 1. 
162 

163 Implementation 

164 Resistance surfaces can be optimized using cost distances (least cost path) or using circuit-based 

165 resistance distances. Cost distances are calculated using the R package gdi stance (van Etten 

166 2014), while electrical current resistances are calculated using the free, open-source software 

167 CIRCUITSCAPE (v4.0 or greater; McRae & Shah 2009). If using ORCUITSCAPE, input data 

168 formats must conform to those of CIRCUITSCAPE (McRae & Shah 2009). Currently, 

169 CIRCUITSCAPE can only be executed from R through Resi stanceGA on computers using 

170 Windows operating systems. Prior to running any of the optimization functions (MS_opti m, 
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171 SS_optim, Resi stance .Opt_multi or Resistance. Opt_single) the 

172 CS . prep/gdi st . prep , and GA. prep preparation functions must be run. These functions create 

173 and format the inputs and data objects necessary to run CTRCUITSCAPE/gdi stance, 

174 parameterize the genetic algorithm, and fit MLPE mixed effects models. These functions have 

175 been developed to require a minimum input from the user. For instance, at minimum, all that 

176 needs to be specified to run GA . prep is the directory where the .asc files to be optimized reside. 

177 However, all available arguments of the ga function can be set by the user to modify the genetic 

178 algorithm. The arguments and default settings of the preparation functions are described in Table 

179 2. 

180 There is no established framework for optimization of resistance surfaces, but an 

181 envisioned work flow is detailed in Figure 2. Genetic algorithms are stochastic optimization 

182 procedures; therefore it is highly advised to run all optimizations at least twice to confirm 

183 convergence and parameter estimates. Also, because boundaries are placed on the parameter 

184 space searched, if the optimized resistance values are at or near the limits set, the optimization 

185 should be rerun after expanding the search space. It is important to note that simultaneous 

186 optimization of multiple surfaces often results in parameter values that differ from truth. 

187 However, the relative resistance values of the resultant resistance surface are highly or perfectly 

188 correlated with truth (see worked example), and the resistance surfaces are identical once they 

189 are rescaled to a minimum resistance of one. Users of Resi stanceGA should be aware of this 

190 fact, and interpret results accordingly. 

191 Genetic algorithms are effective at finding an optimal solution, but the can be 

192 computationally intensive. To ensure that parameter space is adequately searched, the population 

193 of individuals produced each generation must be of sufficient size. In Resi stanceGA the default 
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194 setting is to produce a population that is 15 times the number of parameters being optimized, up 

195 to a maximum population size of 100 individuals. For example, when a continuous surface is 

196 being optimized, 45 individuals (3 parameters x 15) will be produced each generation, and a 

197 typical optimization takes 50-300 generations. This results in the creation of 2 250-13 500 

198 resistance surfaces and CIRCUITSCAPE/gdi stance runs. Therefore, the greatest impediment 

199 to optimizing resistance surfaces is time. Both the spatial extent and the number of sample 

200 locations in the analysis affect the time necessary to calculate pairwise resistance or cost 

201 distances. On a computer with an Intel i7 3.4 GHz processor, 25 sample locations distributed 

202 across a 250 x 250 cell resistance surface takes ~5 seconds to complete an iteration with 

203 CTRCUITSCAPE. The same surface with 100 sample locations takes -14 seconds. Finally, 25 

204 sample locations on and a 1 000 x 1 000 cell surface can take up to 90 seconds to complete. 

205 Because CIRCUITSCAPE calculates resistance over all possible pathways, it is more 

206 computationally intensive than calculating least cost paths. On average, optimization with 

207 gdi stance is about three times faster than calculating resistance distances with 

208 CIRCUITSCAPE. To further reduce the optimization time when using least cost paths, the 

209 optimization can be run in parallel by setting paral 1 el = TRUE in GA. prep . Additional 

210 strategies to reduce the runtime of CIRCUITSCAPE/gdi stance include modifying the 

211 connection scheme (Neighbor. Connect/di recti ons) from the default setting of 8 to 4, and 

212 reducing the resolution of the resistance surfaces (i.e. increase cell size). The latter will reduce 

213 the total number of cells on the surface, and generally has limited effect on the quantification of 

214 relative resistances across the landscape (McRae et al. 2008). 
215 

216 Worked examples 
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217 The examples below use simulated data provided with Resi stanceGA. Code to generate the 

218 example resistance surfaces and spatial point data can be found in Appendix SI. More extensive 

219 examples and further descriptions of the functions included in Resi stanceGA can be found in 

220 the vignette included with the package. 
221 

222 Example 1: Single surface optimization 

223 This example optimizes a single resistance surface using CTRCUITSCAPE. 

224 # Get data 

225 # install 'devtools' package, if needed 

226 if (! ("devtools" %\ n% 1 i st . f i 1 es( . 1 i bPathsO))) { 

227 install . packages("devtool s") 

228 } 
229 

230 1 i brary(devtool s) 

231 i nstal l_gi thub("wpeterman/Resi stanceGA") # Download package 

232 1 i brary(Resi stanceGA) # Load package 

233 rm(list = ls()) 
234 

235 # Create directory for reading/writing CIRCUITSCAPE files and results 

236 if ("ResistanceGA_Examples"%in%di r("C:/")==FALSE) 

237 di r. c reate (file. path ("C:/", "Resi stanceGA_Examples")) 
238 

239 # Create a subdirectory for the first example 

240 di r .c reate (fi le . path("C:/ Resi stanceGA_Exampl es/" , "Si ngl eSurface")) 
241 

242 # Directory to write .asc files and results 

243 write. dir <- "C:/ Resi stanceGA_Exampl es/Si ngl eSurface/" 
244 

245 # Load example data, write continuous surface as \asc~ file 

246 data(resi stance_surfaces) 
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247 continuous <- resi stance_surfaces [ [2] ] 

248 wri te Raster (conti nuous , f i 1 ename=pasteO(wri te . di r , "cont . asc") , overwri te=TRUE) 
249 

250 # Load sample location data 

251 data(sampl es) 

252 wri te . tabl e(samples , 

253 f i 1 e=pasteO(wri te . di r , "sampl es . txt") , 

254 sep="\t" ,col . names=F, row.names=F) 
255 

256 # Create a spatial points object for plotting 

257 sampl e . local es <- Spati al Poi nts(sampl es [ , c(2 , 3)] ) 
258 

259 # Run preparation functions 

260 GA. inputs <- GA.prep(ASCII .di r=write.di r , 

261 max.cat=500, 

262 max.cont=500, 

263 seed=555, 

264 quiet=TRUE) 
265 

266 CS. inputs <- CS . prep(n . POPS=l ength(sampl e . 1 ocal es) , 

267 CS_Poi nt . Fi 1 e=pasteO(wri te . di r , "sampl es . txt") , 

268 CS . program= 1 "C : /Prog ram Fi 1 es/Ci rcui tscape/cs_run . exe" 1 ) 
269 

270 # Transform the continuous surface 

271 r.tran <- Resi stance . tran(transformation="Monomol ecular" , 

272 shape=2, 

273 max=275, 

274 r=conti nuous) 
275 

276 # Visualize the transformation (Fig. 3) 

277 plot.t <- Plot. trans(PARM=c(2, 275) , 

278 Resi stance=conti nuous , 
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279 transformati on="Monomol ecul ar") 
280 

281 # Calculate pai rwi se resistance values to use as the response 

282 CS. response <- Run_CS(CS . i nputs=CS . i nputs , 

283 GA. i nputs=GA. i nputs , 

284 r=r.tran) 
285 

286 # Rerun CS.prep" to include response 

287 CS. inputs <- CS . prep(n . POPS=l ength(sampl e . 1 ocal es) , 

288 response=CS . response , 

289 CS_Point . Fi 1 e=pasteO(wri te .di r , "sampl es . txt") , 

290 CS . program=' "C:/Program Fi 1 es/Ci rcui tscape/cs_run . exe" 1 ) 
291 

292 # Run optimization 

293 SS_RESULTS <- SS_opti m(CS . i nputs=CS . i nputs , 

294 GA.inputs=GA. inputs) 
295 

296 # Compare results 

297 SS_table <- data.f rame(c("Monomolecular" , 2.0, 275), 

298 t(SS_RESULTS$Conti nuousResul ts [c (3 : 5) ] ) ) 

299 colnames(SS_table) <- c("Truth", "Optimized") 

300 SS_table 

301 Truth Optimized 

302 Equation Monomol ecul ar Monomol ecul ar 

303 shape 2 1.999999 

304 max 275 274.9982 
305 

306 The exact transformation parameters were recovered in 147 iterations (Fig 3). 

307 Example 2: Multisurface optimization 

308 This example simultaneously optimizes three resistance surface using least cost paths. 

309 # Get data 
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310 data(resi stance_surfaces) 

311 data(sampl es) 
312 

313 # Create a spatial points object 

314 sampl e . local es <- Spati al Poi nts(sampl es [ , c(2 , 3)] ) 
315 

316 # Create a subdirectory for the second example 

317 di r .create(fi le . path("C:/Resi stanceGA/" , "MultipleSurfaces")) 
318 

319 # Run ~GA.prep~ 

320 GA. inputs <- GA . prep (ASCII . di r=r esi stance_sur faces , 

321 Results. di r="C:/ResistanceGA/MultipleSurfaces/" , 

322 max.cat=500, 

323 max.cont=500, 

324 seed = 321, 

325 parallel = 4) # Run on i n parallel on 4 cores 
326 

327 # Run ~CS.prep~ functions 

328 gdi st . i nputs<-gdi st . prep(n . POPS=l ength (sample . local es) , 

329 sampl es=sample. locales) 
330 

331 # Set parameters for transforming and combining resistance surfaces 

332 PARM <- c(l, 250, 75, 6, 3. 5, 150, 1,350) 
333 

334 # parm<- c(l, # First feature of categorical 

335 # 250, # Second feature of categorical 

336 # 75, # Third feature of categorical 

337 # 6, # Transformation equation for continuous surface 

338 # 3.5, # Shape parameter 

339 # 150, # Scale parameter 

340 # 1, # First feature of feature surface 

341 # 350) # Second feature of feature surface 
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342 

343 # Combine resistance surfaces 

344 Resist <- Combi ne_Surfaces(PARM=PARM, 

345 gdi st . i nputs=gdi st . i nputs , 

346 GA. i nputs=GA. i nputs , 

347 OUt=NULL, 

348 rescale=TRUE) 
349 

350 # Create the response surface 

351 gd. response <- Run_gdi stance(gdi st . i nputs=gdi st . i nputs , 

352 GA. i nputs=GA. i nputs , 

353 r=Resist) 
354 

355 # Re-run ~gdist.prep~ to include the response 

356 gdi st . i nputs<-gdi st . prep(n . POPS=l ength (sample . local es) , 

357 response=l owe r (as . matrix(gd . response)) , 

358 sampl es=sampl e . locales) 
359 

360 # Run ~MS_optirrf 

361 Multi .Surface_optim.gd <- MS_optim(gdi st . inputs=gdi st . i nputs , 

362 GA. i nputs=GA. i nputs) 
363 

364 # view optimized parameter values with the true simulation values 

365 Summary . tabl e <- data.f rame(PARM, round(t(Multi . Surface_optim.gd@solution) ,2)) 

366 colnames(Summary.table)<-c("Truth" , "Optimized") 

367 row. names(Summary . tabl e)<-c("Categoryl" , "Category2", "Category3" , 

368 "Transformation", "Shape", "Max", "Featurel", "Feature2") 

369 Summary . tabl e 

370 Truth Optimized 

371 Categoryl 1.0 1.00 

372 Category2 250.0 268.27 

373 Category3 75.0 80.42 
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379 

380 # Compare the true and optimized resistance surfaces 

381 optim. resist <- Combi ne_Surfaces(PARM=Mul ti .Surface_optim.gd@solution, 

382 gdist. inputs = gdi st . i nputs , 

383 GA. inputs = GA. inputs) 
384 

385 resi st . stack <- stack(Resi st , optim. resi st) 

386 names(resi st. stack) <- c("Truth", "Optimized") 

387 pai rs(resi st . stack) # Correlation plot (Fig. 4) 

388 plot(resi st . stack) # Resistance maps (Fig. 4) 
389 

390 This example took 210 iterations to optimize. As described above, the optimization has 

391 recovered different parameter values than those that were simulated, but the relative values of the 

392 simulated and optimized surfaces are identical (Fig. 4). 
393 

394 Obtaining Resi stan ceGA 

395 Resi stanceGA is hosted on GitHub, and can be downloaded using the i nstal 1 . github 

396 function from the devtool s package: 

397 install . gi thub("wpeterman/Resi stanceGA") 

398 Following download, the package can be loaded: 

399 1 i brary(Resi stanceGA) 

400 View the vignette with more demonstrations of Resi stanceGA functions: 

401 vignette( ' Resi stanceGA' ) 
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402 Alternatively, the HTML vignette can be viewed here: http://goo.gl/n7VOZx 
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Table 1. Functions of the Resi stanceGA package 



Function 



Combi ne_Su rf aces 
CS. prep 

Diagnostic. Plots 

GA. prep 

gdi st . prep 
Grid. Search 

lower 
MLPE. 1mm 
MS_optim 

Plot .trans 

Resi stance .Opt_multi 



Returned objects 



Description 



Resistance. tran 



Run_CS 



R raster object 



Named list 



.tif file 



Named list 



Named list 

contour plot, R data 
object 

R data object 

lmer object 

GA object, diagnostic 
plots, .txt summary 
file 

ggplot object 
AICc value 



Resi stance. Opt_single AJCc value 



R raster object, .asc 
file 



R raster object or R 
data object 



Combine multiple resistance surfaces into 
a new composite surface based on 
specified parameters 
This function prepares and bundles the 
inputs necessary to run CTRCUITSCAPE 

This function returns a multipanel figure 
including the histogram of residuals and 
qqplot from a fitted mixed effects model 

This function prepares and compiles the 
objects and commands needed to execute 
the genetic algorithm 
This function prepares and bundles the 
inputs necessary to run gdi stance 

For a single continuous surface, this 
function can be used to visualize the AICc 
response surface 

Convenience function to obtain the lower 
half of a square distance matrix 

Fits a maximum likelihood population 
effects mixed effects model in 1 me4 

This is a wrapper function for 
simultaneously optimizing multiple 
surfaces. 

Implements and plots a continuous surface 
transformation 

This is the function passed to g a to 
optimize multiple resistance surfaces 
simultaneously 

This is the function passed to g a to 
optimize resistance surfaces in isolation 

Applies specified transformation to 
continuous resistance surface and returns 
a raster object. Optionally, a .asc file can 
be exported 

This function executes CTRCUITSCAPE 
from R and returns either the lower half 
of the pairwise resistance distance matrix 
or the cumulative current map produced 
by CTRCUITSCAPE 
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Run_gdi stance costDistance matrix This function executes gdi stance and 

object from returns the costDistance matrix object 

gdi stance 

SS_opti m .tif files, .csv This is a wrapper function for optimizing 

summary tables surfaces in isolation. All surfaces in a 

common directory will be optimized in 
turn, and numerous summary tables of 
optimized parameters and AICc values are 
produced 
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483 Table 2. Arguments of the required preparation functions and their default settings. Only one of 

484 either CS . prep or gdi st . prep needs to be run, depending upon whether optimization will use 

485 CIRCUITSCAPE or gdi stance. 



Function 


Arguments 


Defaults 


Comments 


CS . prep 


n . POPS 


Must be defined 






response 


NULL 


Must be defined to 
run optimization 




CS_Poi nt . Fi 1 e 


Must be defined 






CS . program 


"C:/Program Files/Circuitscape/cs_run.exe"' 


Default Windows 
installation location 




Neighbor .Connect 


8 




gdi st . prep 


n . POPS 


Must be defined 






response 


NULL 


Must be defined to 
run optimization 




samples 


Must be defined 






transition Function 


function(x) l/mean(x) 






di recti ons 


8 






long 1 al 


FALSE 




GA.prep 


AbLll . en r 


Must be defined 


.asc files should be in 




Mi n . Max 




lllcll uwil UlicLLUiy 




"max" 


Must be "max" when 
optimization with ga 




mi n . cat 


0.0001 






max. cat 


2500 






max. cont 


2500 






cont . shape 


NULL 






pop. mult 


15 






percent . el ite 


0.05 






type 


"real-valued" 


Must be "real-valued" 




pcrossover 


0.85 






pmutation 


0.1 






maxi ter 


1000 






run 


25 






keepBest 


TRUE 






population 


" gareal_Population " 






selection 


"gareal_lsSelection" 






crossover 


"gareal_blxCrossover" 
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mutation "gareal_raMutation" 
seed NULL 
quiet FALSE 



486 
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487 Figure 1. There are eight continuous resistance surface data transformations implemented in 

488 Resi stanceGA . Prior to transformation, the original continuous resistance surface had values 

489 ranging from 0-10. The shape and magnitude of each transformation are each controlled by a 

490 single parameter. All transformations in the figure have a shape parameter value of 3, and 

491 maximum value parameter of 100. Linear relationships are not explicitly incorporated, but all 

492 monomolecular functions become linear as the shape parameter increases. 

493 

494 Figure 2. Flow chart depicting a potential work flow for optimizing resistance surfaces. Analysis 

495 begins by optimizing each resistance surface in isolation. If multiple surfaces are supported (e.g., 

496 AAICc < 4), these surfaces can be simultaneously optimized to create a composite surface. If 

497 desired, different combinations of the best supported single surfaces can be optimized and 

498 compared. Inputs into the optimization are shown in yellow ovals, optimization steps are blue 

499 rectangles, decision points are white diamonds, and final products from optimization are green 

500 ovals. 

501 Figure 3. Optimization of a continuous surface that was transformed using a monomolecular 

502 function following Example 1 in the main text (transformation = Monomolecular, shape = 2.0, 

503 maximum = 275). Optimization with CTRCUITSCAPE took 147 iterations of the genetic 

504 algorithm and precisely recovered the transformation parameters resulting in perfectly correlated 

505 resistance surfaces. 

506 

507 Figure 4. Simultaneous optimization of three resistance surfaces (Example 2 in the main text). 

508 The original categorical, continuous, and feature resistance surfaces are shown on the left and 

509 their actual resistances following transformation/rescaling are shown in the middle. When these 

510 surfaces are added together to create a composite, the minimum resistance value is no longer 1, 

511 so the surface is rescaled to present the relative resistance values (scaled composite resistance). 

512 Optimization using least cost paths took 215 iterations of the genetic algorithm when run in 

513 parallel on 4 cores. The relative resistance values of the optimized surface and the true resistance 

514 surface are perfectly correlated. 
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